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Abstract 

The Dorokhov-Mello-Pereyra-Kumar (DMPK) equation, which describes the 
distribution of transmission eigenvalues of multichannel disordered conduc- 
tors, has been enormously successful in describing a variety of detailed trans- 
port properties of mesoscopic wires. However, it is limited to the regime of 
quasi one dimension only. We derive a one parameter generalization of the 
DMPK equation, which should broaden the scope of the equation beyond the 
limit of quasi one dimension. 
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Quantum transport in a disordered N-channel mesoscopic conductor can be described in 
the scattering approach, initiated by Landauer Jl[], in terms of the joint probability distri- 
bution of the transfer matrices Under very general conditions based on the symmetry 
properties of the transfer matrices and within the random matrix theory framework [|| , the 
joint probability density of the transmission eigenvalues can be expressed as an evolution 
with increasing length of the system according to a Fokker-Planck equation known as the 
Dorokhov-Mello-Pereyra-Kumar (DMPK) equation PJJ51. Such a random matrix approach 
has been found to be very useful in our understanding of the universal properties in a wide 
variety of physical systems in condensed matter as well as nuclear and particle physics ||. 
In particular, the DMPK equation has been shown to be equivalent |7|,[J to the descrip- 
tion of a disordered conductor in terms of a non-linear sigma model |§ obtained from the 
microscopic tight binding Anderson Hamiltonian for non interacting electrons, and is con- 
sistent with perturbative calculations and experiments |T0|j3|JTT|| . The equation has been 
solved exactly |T!2fl , and level correlation functions can be obtained J7| using the method of 
biorthogonal functions [[HJ. Because it is extremely difficult to evaluate any higher order 



correlation function in the sigma model approach, the DMPK equation is more suitable 
to study the conductance distribution in mesoscopic systems. In recent years it has been 
applied to a variety of physical phenomena, including conductance fluctuations, weak local- 



ization, Coulomb blockade, sub-Poissonian shot noise, etc. [llj. One major limitation of the 



DMPK equation however is that it is valid only in the regime of quasi one dimension (ID), 



where the length of the system is much larger than its width fl^ , [TT]| . While the dependence 
on geometry of some of the transport properties have been obtained perturbatively Jl5| in 
the metallic regime, only limited progress have been made on the extension of the DMPK 



equation to higher dimensions [16,17]. Currently, there exists no theory for the statistics 
of transmission levels for all strengths of disorder beyond quasi ID. This is a particularly 
severe shortcoming; the important question of the nature of the expected novel kind of uni- 
versality of the distribution of conductance near the metal-insulator transition ]T8|jr9| can 
not be studied within the powerful DMPK framework, because the transition exists only in 
higher dimensions. 

In this work we argue that the generalization of the DMPK equation to higher dimen- 
sions require the relaxation of certain approximations made in the derivation, and suggest 
a phenomenological way to implement them within the random matrix framework. This 
allows us to obtain a simple generalization using a phenomenological parameter and the 
conservation of the total probability. We obtain corrections to the mean and variance of 
conductance as a function of the parameter using the generalized equation and discuss the 
implications of the results. We argue that the generalized equation should be valid beyond 
quasi one dimension. 

In the scattering approach, the conductor of length L is placed between two perfect leads 
of finite width. The scattering states at the Fermi energy define N channels. The 2N x 2N 
transfer matrix M relates the flux amplitudes on the right of the system to that on the left 
PJ3[]. Flux conservation and time reversal symmetry (in this paper, for simplicity, we will 
restrict ourselves to the case of unbroken time reversal symmetry only) restricts the number 
of independent parameters of M to N(2N + 1), and can be represented as M 
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where u, v are N x N unitary matrices, and A is a diagonal matrix, with positive elements 
Xi,i = 1,2...N. The physically observable conductance of the system is given by g = 
5^(1 + Aj)" 1 . Thus the distribution of conductance can be obtained from the distribution 
of the variables A,. 

In order to understand the nature of the approximation used in DMPK and to motivate 
our generalization, we will briefly review the derivation of DMPK following ref. Q. In 
this approach, an ensemble of random conductors of macroscopic length L ^> I, where I is 
the mean free path, is described by an ensemble of M random matrices, whose differential 
probability depends parametrically on L and can be written as dP L (M) = p L (M)dp,(M). 
Here dfi(M) is the invariant Haar measure of the group, given in terms of the parameters 
in (1) by 



d(i(M) = J(A) 



N 



lid\ 



d[i(u)dfi(v) 



(2) 



where 



■/(A) = n>- A ;l ( 3 ) 

i<j 

and dfi(u) and dfi(v) are the invariant measures of the unitary group U(N). When a con- 
ductor of length Lq described by a transfer matrix Mo is added to a conductor of length 
Li and transfer matrix Mi to form a conductor of length L = L\ + L and transfer matrix 
M = Mo Mi, the probability density pi(M) satisfies the combination rule 

(Pl 1+ l (M)) Lo = J p Ll (A^ - 1 K (M )rf/i(M ), (4) 

where the angular bracket represensts an ensemble average. For Lq <C /, the small change 
in the transfer matrix can be expected to lead to a small change in the parameters Aj, and 
one can expand the probability density as 

(pl 1+ l (X)) Lo = (Pii(A + SX)) Lo = PLl (X) + E (5X a ) Lo + \ £ (SX a 5X b ) Lo 

(5) 

Since the changes in A a are small, we can use perturbation theory to evaluate their averages. 
We can also expand the left hand side in powers of Lq. The resulting equation, keeping only 
terms first order in Lq on the left hand side, is given by 

L o^ = E(i + 2 ^)^ (E Kv' c ;< a ) L + E Vi + A a )|| (^2 K0- + K)WS) L 

A a + Xb + 2A a A& dp L-^ , , I*,*,, \ , . 

+ E x v JTT-il^K^ + KKa V cb V cb V ca ) . (6) 

a+b A « - A b OA a \ c I Lo 

Here the primed variables correspond to the added small conductor of length L . 
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The above equation (6) is quite general. It is based on the symmetry properties of the 
transfer matrices and the combination principle for adding two conductors. These principles 
should remain valid beyond quasi one dimension. It is the further approximations on the 
averages in equation (6) made in deriving DMPK that limits DMPK to quasi one dimension. 
There are two major approximations involved: 

(i) The 'isotropy' assumption is used to decouple the averages over the products of the 
parameters A and the unitary matrices v. Once decoupled, the averages over the products 
of the unitary matrices alone can be explicitly obtained to give 

/ * / \ J_. //*/*/ / \ 1 . /i / |4\ ^ /y\ 

V ca V ca/- N 'i \ V ca V cb V cb V ca/ ~ + 1 y 1 / ~ N + 1 ) ' ^' 

while the average over the trace of \' c is taken to be proportional to L . In particular, 
(J2cK) = NL /l, where / is the mean free path, consistent with the Born approximation 
for the transmission amplitude valid for small L . 

(ii) The second approximation is based on the expectation that the averages of the 
products of A^, are higher orders in L , and therefore negligible. In particular, this means 
that the terms proportional to J2 C K 2 are neglected in equation (6). 

The above two approximations, together with the identity 

lead to the well known DMPK equation: 

d P (\y 



dp _2 1 v d 

d(L/l) ~ N + l J(A) \ W a 



\ a (l + X a )J(\y 



a 



(9) 



where J(A) is defined in (3). 

We will first show that beyond quasi one dimension, the second approximation fails, 
namely J2 C K 2 ^ s °f ^ ne same order in L as J2 C K an< i therefore can not be neglected. In 
this case we will show that the total probability can not be conserved within the decoupling 
approximation. We will then introduce phenomenological parameters for the averages over 
the products in (6), and show that the conservation of total probability require a very specific 
generalization of the DMPK equation involving a single additional parameter. Finally we will 
evaluate the corrections to the mean and variance of the conductance using the generalized 
DMPK as a function of the parameter and interpret the results. 

To go beyond quasi ID, we start with a conductor of length Lq along x and width W 
along y and z, with scattering potential V(x,y,z). To see how the second approximation 
fails, we will consider, for simplicity, a square well potential adequately approximated by 
a repulsive delta function at x — 0, i.e V(x,y,z) = Vr(y, z)S(x). Writing the Schrodinger 
wavefunction as ty(x,y,z) = J2i i>i(y, z )(pi( x ), where ipi(y,z) are the transverse eigenfunc- 
tions in the perfectly conducting lead, chosen to be real, we obtain the system of coupled 
equations for the N channels 

<f>'!{x) + k 2 Mx) = E Ka{x)<l>j{x), (10) 
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where the prime denotes a derivative with respect to X . k'l £1X6 the wavevectors in channel i, 
and Kij are the coupling constants given by Kij(x) = (2m /h) J J dydzipj(y, z)Vr(y, z)ipi(y, z). 
We are interested in the transfer matrix M that connects the solution on the left side of the 
conductor with that on the right side. The transfer matrix satisfying the flux conservation 
and time reversal symmetry can be written in the form 

«=( 1+ S * ,), (id 



A* 1 + A* 

where 1 and A are N x N matrices and A^- = Kij/2iki. Note that A is pure imaginary 
but not symmetric. The parameters A that satisfy the DMPK equation in quasi ID are 
the eigenvalues of the matrix X = [Q + Q^ 1 — 2 ■ l]/4, where Q = M^M ||. From flux 
conservation, Q^ 1 = S Z QS Z where S 2 is the third Pauli matrix with 1 and replaced by 
(N x N) 1 and matrices. It is easy to see that X is block diagonal, each block given by 
a sum of two matrices Xi = (A + A^)/2 and X 2 = A^A. The important point is that Xi 
is traceless, so tr(Aj) is given by tr(X 2 )=tr(AtA). On the other hand, X\ does contribute 
to tr(Af)=tr(Xx + X 2 ) 2 , where tr(Xx) 2 =tr(A 2 + At 2 + A+A + AAt)/4. Clearly it is of the 
same order as tr(Aj), and can not be neglected. 

It is now straightforward to show that keeping the tr(A 2 ) terms in (6) and using the 
decoupling approximation of the averages of v and A lead to a breakdown of the conservation 
of total probability. Suppose (j2 c K 2 ) = a Lo/l- Then using (7) for the averages over t>, we 
get a correction term to the DMPK equation equal to 

Clearly this is not a sum of total derivatives and the resulting equation does not conserve 
total probability |20 . 



It is therefore clear that in order to go beyond quasi ID, we need to relax both ap- 
proximations. We propose a simple phenomenological way to take care of both. Instead 
of computing the three averages in (6) explicitly, we start with the following very general 
ansatz: 

(y2K V 'ca* V ca) = yl (J2K( 1 + K) V 'ca V 'cb* V cb V 'ca) = YW+l^ 

where fii and fi 2 are arbitrary dimensionless parameters, which can be functions of N. 
Clearly, \i\ = /i 2 = 1 gives back the quasi ID limit. Note that any additional parameter in 
the first term will only serve to redefine the mean free path, so there are only two additional 
parameters possible. With this ansatz, equation (6) becomes 

dp N — 1 „ 9p 2/x 2 /1 d 2 p 



d(L/l) ^ ^JV + l'V^ ' "" a 'd\ a ' iV + l^" a ^ ' " a 'd\l 

2ui , , , , 1 dJ dp . <N 

+wh^^iwA (14) 
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We now demand that the parameters /xi and n% are such that the right hand side can be 
written as a sum of total derivatives in order to ensure the conservation of total probability. 
Note that the special choice fix = ^ — 1 makes the coefficients of all three terms on the 
right hand side of (14) the same, and then the three terms can be written as a sum of 
derivatives after multiplying by J(X). It may appear at first that with two parameters and 
three terms, no other choice is possible, except for a trivial multiplicative factor for all three 
terms which can be absorbed in the redefinition of the mean free path. However, we note 
that if we choose 



N -1, 



N + 1 J iV + 1' 



together with a renormalization of the measure 



J — > J 7 ; 



01 

7 = — , 

A*2 



(15) 



(16) 



then (14) can be rewritten as 
dp 2 



d{L/V) 



N+lJ^X)^dX a 



E 







X a {l + X a )J\X) 



dp{\) 

dX n 



(17) 



where Z' = Z//i2 is a renormalized mean free path. Equation (17) is our one parameter gener- 
alization of the DMPK equation (9), where the parameter 7 enters in the renormalization of 
the measure as in (16). Note that in the absence of time reversal symmetry or in the presence 
of spin-orbit scattering, the measure is changed in a similar way by an exponent f3 = 2, 4 
respectively |14j2l|. However, in our present case with time reversal symmetry, (3 = 1, and 
the exponent 7 is in general non integral. Clearly 7 = 1 is the quasi ID limit. From the 
relation between fii and fi2, and the condition that both fix and \i2 must be positive, we find 
the following restrictions: 



N + l n N+l 
0< ^< < 02 < — ^— 



(18) 



This means that the only restriction on the parameter 7 is that it is positive. In general, it 
can be a function of N. 

We can try to interpret the phenomelogical parameter 7 by comparing with known 
results. The expectation value of any function F(\), defined as 



(F) 



{L/V) 



N 



F(x)p {L/n (x)r(\)i[dx a: 



(19) 



a=l 



follows an evolution equation which can be obtained by multiplying both sides of (17) by 
J 7 (A)F(A) and integrating over all X a , giving 



9{F) t 
ds 



E 



dF d 2 F 
l + 2X a )— + X a (l + X a ) w 



dX n 



2 a+b X a ~ X b 



dF 



(20) 
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where s = L/V . If 7 is independent of N, then we can use the method of moments in 



14] to obtain the average and variance of the conductance g = + Aj) 1 as a power 



series in s/N <C 1 in the large N and large s limit. We find that to leading order, (g) = 
Nl'/L — (2 — j)/3j, and var(g) = (g 2 ) — (g) 2 = 2/157 . As expected, 7 = 1 gives back the 
quasi ID results. However, in general the variance decreases with increasing 7. Comparing 



with the result var(g) ~ ^jL y L z /L for a rectangular conductor with length L x = L and 
cross-section L y L z [p2fl , we see that the parameter 7 can be identified with the aspect ratio 



Jy L z in this diffusive transport regime. 
If 7 ~ 1/-/V, obtained for \xi = za/V, which is consistent with the restriction (18) if 
v < 1/2, then we need to divide both sides of (20) by 7, so that the renormalized mean 
free path I" = I'/j = l/fii is independent of N. Then assuming v <C 1, it is possible to 
obtain corrections to the 1/N expansion up to linear order in v using the above method of 
moments. The result is that the corrections are larger by a factor vfiis, which signals the 
breakdown of the expansion in the large s limit. It would be interesting to obtain a more 
rigorous solution of (17) for arbitrary 7. 

In summary, by relaxing certain approximations in the derivation of the DMPK equation 
(9) which limits it to the quasi ID regime only, we have derived a one parameter general- 
ization given in eq. (17), based on a phenomeno logical ansatz and the conservation of total 
probability. The geometry dependence of the parameter, obtained in the diffusive limit by 
evaluating the correction to the variance of the conductance beyond its quasi ID value, sug- 
gests that the generalized equation should be applicable beyond the quasi ID regime. This 
should broaden the scope of the DMPK approach. 
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